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ABSTRACT 


We present arcminute-resolution intensity and polarization maps of the Galactic center made with the 
Atacama Cosmology Telescope. The maps cover a 32 deg? field at 98, 150, and 224 GHz with |I| < 4°, 
|b| € 2°. We combine these data with Planck observations at similar frequencies to create coadded 
maps with increased sensitivity at large angular scales. With the coadded maps, we are able to resolve 
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many known features of the Central Molecular Zone (CMZ) in both total intensity and polarization. 
We map the orientation of the plane-of-sky component of the Galactic magnetic field inferred from 
the polarization angle in the CMZ, finding significant changes in morphology in the three frequency 
bands as the underlying dominant emission mechanism changes from synchrotron to dust emission. 
Selected Galactic center sources, including Sgr A*, the Brick molecular cloud (G0.2534-0.016), the 
Mouse pulsar wind nebula (G359.23-0.82), and the Tornado supernova remnant candidate (G357.7- 
0.1), are examined in detail. These data illustrate the potential for leveraging ground-based cosmic 
microwave background polarization experiments for Galactic science. 


1. INTRODUCTION 


Some of the most extreme interstellar environments in 
the Galaxy are found in the Galactic center (e.g., Bat- 
tersby et al. 2020). The inner ~ 500 pc of the Milky Way 
is home to the Central Molecular Zone (CMZ), the dens- 
est concentration of molecular gas in the Galaxy, with 
a mean density of ~ 10* cm-? (Güsten 1989; Ferrière 
et al. 2007). The surface density of dense gas greatly 
exceeds that found in nearby star-forming molecular 
clouds, with the average gas surface density transition- 
ing from ~ 5Mg pc”? to several hundreds My pc”? as 
one reaches the inner 200 pc of the Galaxy (see Morris 
& Serabyn 1996, for a review). Standard prescriptions 
predict that the CMZ should be an extremely active site 
of star formation, and yet the observed star formation 
rate is low; by some estimates, an order of magnitude 
or more below predictions (e.g., Longmore et al. 2013; 
Barnes et al. 2017; Nguyen et al. 2021, and references 
therein). 

'The apparently inefficient star formation in the CMZ 
makes this region an ideal testbed for star formation 
theories, with many factors proposed to explain the ob- 
servations. These include the strong magnetic field in 
the Galactic center (Crutcher et al. 1996; Chuss et al. 
2003; Ferriére 2011), the rate of mass inflow to the CMZ 
(Sormani & Barnes 2019), the strength and compress- 
ibility of turbulence in the CMZ (Federrath et al. 2017), 
and the possibility that we are observing a relatively 
quiescent period between episodic bursts of star forma- 
tion (Kruijssen et al. 2014; Krumholz & Kruijssen 2015). 
Furthermore, the CMZ is in some respects a nearby ana- 
log of nuclear rings in other galaxies, including high- 
redshift starbursts. The Galactic center is thus an op- 
portunity for up-close study of the physics relevant to 
the cosmic history of star formation (Kruijssen & Long- 
more 2013; Ginsburg et al. 2019). 

The magnetic field in the vicinity of the Galactic 
center has long been studied with radio polarimetry 
(Ferriére 2009; Morris 2015). The so-called nonthermal 
radio filaments — thin strands of radio-frequency emis- 
sion — were some of the earliest observations to shed light 
on the magnetic field structure toward the Galactic cen- 


ter. The nonthermal radio filaments are, for the most 
part, strikingly perpendicular to the Galactic plane, and 
the intrinsic magnetic field inferred from the Faraday 
de-rotated polarized synchrotron emission tends to lie 
parallel to the long axis of these filaments (Morris & 
Yusef-Zadeh 1985; Yusef-Zadeh & Morris 1987a; Lang 
et al. 1999). 

Polarized dust emission provides a complementary 
means of probing the magnetic field structure in the 
CMZ. Interstellar dust grains emit partially polarized 
thermal radiation because they are aspherical and pref- 
erentially align their short axes parallel to the ambient 
magnetic field (Purcell 1975). The polarization angle 
of the dust emission is thus a line-of-sight (LOS) inte- 
grated probe of the plane-of-sky component of the mag- 
netic field orientation. Polarized dust emission has been 
measured at high angular resolution in small regions to- 
ward a number of CMZ molecular clouds (e.g., Novak 
et al. 2000, 2003; Chuss et al. 2003; Matthews et al. 
2009; Roche et al. 2018). Recently, the balloon-borne 
experiment PILOT presented a 240 um map of the po- 
larized dust emission over the entire CMZ region at 2.2' 
resolution (Mangilli et al. 2019), along with comparisons 
to the lower-resolution 353 GHz polarization data mea- 
sured by the Planck satellite (Planck Collaboration Int. 
XIX 2015). 

The Atacama Cosmology Telescope (ACT) measures 
the polarized microwave sky with higher angular reso- 
lution than the Planck satellite and greater sensitivity 
on small scales. In this paper we present new dedicated 
maps of the Galactic center in total intensity and linear 
polarization in three ACT frequency bands. We combine 
the ACT data with Planck data to augment the map sen- 
sitivity on larger angular scales. The frequency coverage 
of the maps presented here probe a range of physical 
emission mechanisms, enabling a comprehensive view of 
the Galactic center environment. In polarization these 
maps probe both polarized dust and synchrotron emis- 
sion, and in total intensity the maps additionally show 
features from free-free emission and molecular line emis- 
sion from transition frequencies that fall within the ACT 
passbands. These data illustrate the potential of sensi- 


tive cosmic microwave background (CMB) polarization 
experiments for Galactic science. 

We describe the ACT observations in Section 2 and 
the mapmaking and Planck coadd procedures in Sec- 
tion 3. In Sections 4 and 5, we present the maps in 
total intensity and polarization, respectively, and dis- 
cuss derived properties, including emission mechanisms, 
magnetic field orientation, and polarization fraction. In 
Section 6, we identify notable Galactic center objects 
and compare to observations at other frequencies. We 
conclude in Section 7. 


2. OBSERVATIONS 


ACT is a 6-meter off-axis Gregorian telescope located 
at an elevation of 5190m on Cerro Toco in the At- 
acama Desert in Chile (Fowler et al. 2007; Thornton 
et al. 2016). ACT scans the millimeter-wave sky with ar- 
cminute resolution, complementary to the full-sky lower 
angular-resolution measurements from satellite missions 
such as the Wilkinson Microwave Anisotropy Probe 
(WMAP; Bennett et al. 2013) and Planck (Planck Col- 
laboration I 2014). 

In 2019, the target ACT observing fields were ex- 
panded to include the Galactic center region. Between 
2019 June 6 and November 29, a total duration of ~ 35 
hr of data were taken with three Advanced ACTPol 
dichroic detector arrays PA4, PA5, and PA6 (Hender- 
son et al. 2016; Ho et al. 2017; Choi et al. 2018), at 
three frequency bands f090, f150, and f220 centered 
roughly at 98, 150, and 224 GHz, respectively. The beam 
full-width-half-maximum (FWHM) at each band is 2.0’, 
1.4’, and 1’, respectively. The observation field extends 
roughly from —89? to —97? in declination and —33? to 
—25° in right ascension. This study focuses specifically 
on a 32 deg? field near the CMZ with Galactic longitude 
||| € 4? and Galactic latitude |b| < 2°. 

In this paper we present the maps made using the 
nighttime observations only, which constitute roughly 
two-thirds of the total data collected. The daytime ob- 
servations are affected by a time-dependent beam defor- 
mation due to the heating from the Sun that is chal- 
lenging to correct for in detailed high-resolution maps, 
and hence those data are excluded from this analysis. 
Correcting for this beam deformation will be a subject 
of future study, and the daytime observations may be 
included for future versions of these maps. 


3. MAPMAKING 
3.1. Mapmaking with ACT 


'The instrument records observations in the form of 
time-ordered data (TOD) in units of ~ 10 minutes. We 
largely follow the mapmaking pipeline as described in 
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Aiola et al. (2020) with a few key differences, as we 
briefly summarize below. 

First, we cut bad samples affected by glitches in each 
TOD. To prevent bright sources in the Galactic center 
region from being mistaken for glitches, we mask sources 
brighter than 5 mK with a radius of 3’ prior to applying 
the glitch finder. The 5 mK flux limit is chosen such 
that it is low enough to prevent bright sources from be- 
ing mistaken as glitches, but high enough to ensure only 
a tiny fraction of sky is masked. We also note that this 
mask is only applied when identifying glitches and not 
used during mapmaking. Timestreams with outlying 
statistical properties in terms of noise levels and opti- 
cal responsiveness are then flagged and removed from 
the analysis. We further split the dataset into two in- 
dependent subsets for each frequency band and detector 
array respectively, resulting in 12 datasets in total. We 
then obtain the sky maps for each dataset by solving the 
mapping equation, 


d= Pm+n, (1) 


for a set of Stokes parameters (I, Q, U), where d is 
the pre-processed time-streamed data, P is the pointing 
matrix, m is the output map of interest, and n is the 
noise model. This equation yields a maximum-likelihood 
solution for m by inverting 


(PT N-!P)m = PTN-1d, (2) 


where N is the detector-detector noise covariance. 
'There are two notable differences between the pipeline 
used in this study and that presented in Aiola et al. 
(2020). First, we have adopted a new calibration 
method that improves gain stability and reduces biases 
from thermal contamination as compared to the method 
in ACT Data Release (DR) 4 (Aiola et al. 2020) (see 
Appendix A for more details). The second difference 
relates to the handling of point sources and extended 
hot regions that are common in the Galactic center re- 
gion but uncommon in CMB fields. Directly applying 
the mapmaking pipeline in ACT DRA leads to stripes 
around the bright sources caused by model errors as ex- 
plained in Næss (2019). This happens for two reasons: 
(1) A pixelated map does not capture the sub-pixel be- 
havior of the sky. These residuals are proportional to the 
gradient of the signal across a pixel and are often frac- 
tionally small. However, if the sky is sufficiently bright, 
such as in the brightest parts of the Galactic center, they 
can still end up being large in absolute terms. Since the 
map m in Equation 1 cannot capture these residuals, the 
model forces them to be interpreted as part of the noise 
n. (2) The correlated noise model used in the mapmaker 
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Figure 1. Comparison between Planck-only maps (left column) and ACT+Planck coadded maps (right column) in total 
intensity. Rows from top to bottom correspond to f090, f150, and f220 respectively. Each map extends from |/| < 2°, |b| < 1? 
and is plotted on a logarithmic color scale from 0.3-30 MJysr”* for f090 and f150, and from from 3-100 MJysr”* for £220. See 
Figure 14 in Appendix B for a corresponding plot with ACT-only maps. 


induces a nonlocal response to the sub-pixel noise, lead- 
ing to biases on the scale of the noise correlation length. 
To avoid this problem, we first identify the regions that 
source the strongest model errors, namely, the brightest 
parts of the Galaxy, and then eliminate model errors in 
these pixels by allocating an extra degree of freedom for 
each sample that hit them, as described in Neess (2019). 

A caveat concerning these maps is that the bright 
parts of the Galaxy were not masked when building the 
noise model N. The noise model estimator assumes that 
the time-ordered data is noise dominated (d = n), and 
uses this to measure the noise covariance directly from 
d. This breaks down when the telescope scans across the 
Galactic center, resulting in an overestimate of the noise 


amplitude especially on smaller scales. This has two con- 
sequences: (1) The data are weighted sub-optimally in 
Equation (2), resulting in slightly higher noise. Since 
the maps are strongly signal dominated, this can be ig- 
nored. (2) Because the noise model is contaminated by 
the same signal it is applied to, there is a small loss in 
signal power in the maps; pixels where noise happens to 
add constructively to the signal have more power in d 
than in pixels where they partially cancel. Since we use 
inverse-variance weighting, the latter are up-weighted 
compared to the former. The size of this effect is lim- 
ited because the problematically bright regions make up 
a small fraction of the total samples used to build N. 
We have not measured the precise size of this effect, but 
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Figure 2. A side-by-side comparison between Planck only (left) and the ACT+Planck coadded (right) for f150 in polarized 


intensity. 
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Figure 3. Polarization maps in Stokes Q (left column) and U (right column) in Galactic coordinates and using the IAU 
polarization convention. Top to bottom are the f090, f150, and f220 maps, respectively. 


estimate it to be S 1% based on experience with other 
high-signal-to-noise (S/N) regions, and hence we expect 
it to have negligible impact on the interpretation of the 
maps in this paper. This deficiency will be rectified in 
the upcoming ACT DR6 maps. 

A final known issue requiring mitigation is 
temperature-to-polarization (T-to-P) leakage. ACT 
typically scans a given region of the sky both during its 
rising and setting. As the Galactic center region is at 


relatively low declination, rising scans and setting scans 
are poorly cross-linked (for more information on ACT 
scan strategy see Stevens et al. (2018)). Furthermore, 
the ACT beam is known to leak T-to-P at the percent 
level. This beam leakage effect averages down effectively 
in the nominal CMB maps, which are well cross-linked, 
but in the Galactic center region the T-to-P leakage is 
apparent at a ~ 1-296 level that contaminates the po- 
larization maps in the bright Galactic plane. To reduce 
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the contamination from beam leakage, we build a 2D 
leakage beam model for each dataset using observations 
of Uranus made in the same observation year (2019), 
and de-project the expected T-to-P leakage from the 
polarization maps in each dataset (see Appendix B for 
more details). 

Following these methods, we produced two-way split 
maps of the Galactic center region at 0.5’ resolution in 
Plate Carreé (CAR) projection for each frequency band 
(f090, f150, £220) and detector array (PA4, PA5, PAG) 
resulting in a total of 12 maps. 


3.2. Coadd with Planck 


The large angular scales in the ACT maps are affected 
by atmospheric noise contamination and complicated 
co-variances at large scales. These modes can be re- 
covered, however, by coadding ACT maps with maps 
from Planck, which dominate the S/N at large scales 
£ < 1000. In particular, we have used a similar algo- 
rithm as presented in ACT DR5 (Naess et al. 2020), 
in combination with the Planck High Frequency Instru- 
ment (HFT) maps processed through the NPIPE pipeline 
(Planck Collaboration Int. LVII 2020), which are two- 
way split maps featuring improved noise level and sys- 
tematic control as compared to the previous Planck data 
releases. 

As the coadding algorithm is presented in detail in 
Naess et al. (2020), we only briefly summarize the steps 
and note differences here. First, we re-project the 
Planck maps and noise models from HEALPix! (Górski 
et al. 2005) projection with Nsiqe = 2048 into the ACT 
Galactic center observation footprint in CAR projection 
with 0.5' pixelization using bi-cubic interpolation. We 
have used the same passbands as in Naess et al. (2020) 
and similarly matched the Planck 100 GHz maps with 
ACT f090 maps, 143 GHz with ACT f150, and 217 GHz 
with ACT f220. 'This process assumes that the ACT 
and Planck passbands are equivalent. We note that this 
introduces additional scale dependence to the effective 
band centers (Naess et al. 2020). 'This is expected to 
have negligible impact on the results presented here but 
is relevant for component-separation analysis, which will 
be the subject of follow-up work. 

We then solve for the maximum-likelihood coadded 
maps using a block-diagonal equation 
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Band Planck Dataset ACT Dataset 'Total 
f090 100 GHz f090 PA5+PA6 6 
f150 143 GHz f150 PA4+PA5+PA6 8 
f220 217 GHz f220 PA4 4 


Table 1. Subsets of maps coadded at each frequency band. 
All input maps are two-way split maps. The column “total” 
shows the total number of maps coadded in each band. For 
example, 6 different maps went into making the f090 coadd 
map, consisting of two splits from ACT PA4, ACT PA5, and 
Planck 100 GHz, respectively. 


where m; refers to each individual map, B; refers to its 
corresponding beam transfer function, and m refers to 
the final coadded map with a desired beam Boyt, which 
is the ACT beam in this case. n refers to the map noise, 
which is assumed to be Gaussian and block diagonal 
across individual maps, i.e., individual maps have in- 
dependent noise realizations. Of the noise models pre- 
sented in Naess et al. (2020), we have adopted the con- 
stant correlation noise model, though the choice makes 
little difference in practice as we are considering only a 
small patch of sky with close to uniform noise levels. We 
invert Equation (3) to find a maximum-likelihood solu- 
tion to the coadded map at f090 and f150, respectively. 
Because the PA4 array had a poor detector yield over 
the course of the observation, maps at f220 are treated 
differently from the other two frequencies. The result- 
ing excess noise in the ACT f220 maps leads to a lack 
of convergence when solving for a coadded map through 
a maximum likelihood approach. Therefore, we instead 
perform a straightforward inverse-variance weighting in 
Fourier space to obtain the coadded map at f220 (see 
Appendix C for more details). 

One caveat in using the Planck HFI maps is that a 
cosmic infrared background (CIB) monopole model was 
deliberately included on a per-frequency basis due to a 
lack of sensitivity to the absolute emission level. We 
therefore subtracted the CIB monopole in each coadded 
map following Table 12 in Planck Collaboration et al. 
(20202). 

This procedure yields a total of three coadded maps 
in both temperature and polarization at f090, £150, and 
f220, as summarized in Table 1. We present a side- 
by-side comparison between Planck maps and our three 
coadded maps in total intensity in Figure 1, and a similar 
comparison for polarized intensity for f150 in Figure 2. 
It is apparent that the addition of ACT data signifi- 
cantly improves the angular resolution of the maps in 
both temperature and polarization. The coadded polar- 
ization maps are presented in Figure 3 in Galactic co- 
ordinates. We use the IAU polarization convention, in 
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Figure 4. Multifrequency view of the Galactic center region in both total intensity (upper panel) and polarized intensity 
(lower panel). Red, green, and blue correspond to f090, f150, and f220, respectively. In the upper panel, the maps are scaled 
logarithmically from 0.2 to 2 MJy sr^ ! for £090, from 0.214 to 2.14 MJy sr^ ! for f150, and from 1.15 to 10.15 MJy sr^ ! for f220. 
The polarization maps shown in the lower panel are first smoothed with a Gaussian kernel (FWHM= 3.5’) and then scaled 
linearly from 0 to 1 MJy sr ^! for f090, to 1.79 MJy sr ! for f150, and to 8.2 MJy sr ! for £220. 


which the polarization angle measures 0? toward Galac- 
tic North and increases counter-clockwise (Hamaker & 
Bregman 1996). The ACT Collaboration has adopted 
the IAU convention for all ACT data products since 
DR4. This is in contrast to the COSMO convention 
(Górski et al. 2005) adopted in, e.g., the Planck data 
releases, that is related to the IAU convention via a sign 
flip of Stokes U, i.e., Ucosmo = —Utau.- 

A detailed discussion of these maps is presented in 
Section 4 for total intensity maps and in Section 5 for 
polarization maps. The final coadded maps have median 
noise levels of 36 uK-arcmin at f090, 33 yK-arcmin at 
f150, and 270 »K-arcmin at f220. 


4. TOTAL INTENSITY MAPS 


Figure 1 shows the total intensity maps for both 
Planck-only and the coadded maps for our three fre- 


quency bands (f090, f150, £220). Many prominent fea- 
tures that were obscured or unresolved in the Planck 
maps become apparent with the addition of ACT data, 
and qualitative changes in map morphology with fre- 
quency are evident. The Galactic Center Radio Arc 
(GCRA), a prominent filament in the Galactic center, 
is visible at both f090 and f150 near the center of the 
coadded maps and to a lesser extent at f220, consistent 
with it being a strong source of synchrotron radiation 
Paré et al. (2019). 

'The ACT frequency coverage probes a variety of emis- 
sion mechanisms, including synchrotron, free-free, ther- 
mal dust, and molecular line emission, at different levels 
in each of the three bands. To better visualize the differ- 
ent structures probed at each frequency band, we com- 
bine the coadded maps from three frequency bands into 
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Figure 5. Upper panel: Known radio sources found in the Galactic center region. The background image shows a zoomed-in 
view of the multifrequency 3-color image presented in the upper panel of Figure 4. Lower panel: Annotations of selected radio 
and dusty sources in the multifrequency polarized intensity image (presented in the lower panel of Figure 4). Note that we used 


a smoothing with FWHM=2’ to make objects more visible. 


a multicolor image shown in the upper panel of Figure 4. 
The red, green, and blue image channels represent the 
f090, £150, and f220 maps, respectively, after appropri- 
ate rescaling. The intensity scaling (as detailed in the 
Figure 4 caption) was chosen to highlight structures in 
different bands and to make feature identification easier. 
An annotated zoom-in of the three-color intensity map 
in Figure 4 is provided in the top panel of Figure 5. 
The coherent structures visible in the different col- 
ors of Figures 4 and 5 arise from spatial variations 
in the relative strengths of the various emission mech- 
anisms. The radio spectrum of supernova remnants 
(SNRs) originates primarily from synchrotron emission 
(Weiler & Sramek 1988), and thus objects like the SNR 
candidate G357.7-0.1 (“the Tornado") (Milne 1970) and 
SNR0.9--0.1 (Helfand & Becker 1987) appear reddish 
yellow in Figure 4. Similarly, prominent radio sources, 


including Sgr A*, the GCRA, and Sgr B1 (see, e.g., 
Yusef-Zadeh et al. 1984; Pedlar et al. 1989; Bally et al. 
1991) are strikingly highlighted in this color, consistent 
with their strong synchrotron emission spectrum. Pul- 
sar wind nebulae (PWN), like the Crab Nebula, also 
emit highly polarized synchrotron emission with a flat 
spectral index (Gaensler dz Slane 2006), in contrast to 
SNRs, which generally emit synchrotron with a slightly 
lower polarization fraction and a steeper spectrum. 
Thermal emission from interstellar dust dominates 
Galactic emission at far-infrared /submillimeter frequen- 
cies. Known molecular cloud complexes like Sgr B2 
(G0.667-0.031; e.g., Scoville et al. 1975), Sgr C 
(G359.429-0.090; e.g., Liszt & Spiker 1995), and dense 
molecular clouds like the Brick (G0.253+0.016; e.g., 
Longmore et al. 2012), the 20kms”* Cloud (G359.889- 
0.093) and 50km s^! Cloud (G0.070-0.035; e.g., Güsten 


& Downes 1980), and the Three Little Pigs (G0.145- 
0.086, G0.106-0.082, and G0.068-0.075; see, e.g., Bat- 
tersby et al. 2020, for an overview of these molecular 
clouds) thus appear bright blue/green in Figure 5. 

In general, however, the presence of strong molecular 
line emission in the CMZ precludes the simple interpre- 
tation that low frequencies correspond to synchrotron 
emission and high frequencies correspond to dust emis- 
sion. Even in the relatively broad Planck and ACT 
passbands, line emission can dominate the total inten- 
sity in the Galactic center maps. Indeed, Planck Col- 
laboration X (2016) found that 88.6 GHz HCN emis- 
sion can alone account for up to 2396 of the total 
intensity in the Planck 100 GHz band in this region. 
CO(1-0) at 115.3 GHz and CO(2-1) at 230.5 GHz con- 
tribute significantly to the observed emission in the 
Planck 100 and 217 GHz bands, respectively (Planck 
Collaboration X 2016), while other lines such as HCO* 
(89.2 GHz), CS (98.0, 147.0, and 244.9 GHz), !?CO(1-0) 
(110.2 GHz), CN (113.2, 113.5 GHz), H3CO (140.8 and 
218.2 GHz), NO (150.2, 150.5 GHz), SiO (217.1 GHz), 
SO (219.9 GHz), and CO(2-1) (220.4 GHz), among 
others, are also known to be present in the Galactic cen- 
ter (e.g., Liszt & Turner 1978; Sandqvist 1989; Kramer 
et al. 1998; Lang et al. 2002; Takekawa et al. 2014; Pound 
& Yusef-Zadeh 2018; Lu et al. 2021; Schuller et al. 2021) 
and will contribute to the observed emission in the ACT 
and Planck frequency channels. 

The very bright CO(1-0) emission poses a particu- 
lar challenge for our analysis, as it falls comfortably 
within the Planck 100 GHz passband but largely out- 
side that of ACT f090 (see Naess et al. 2020, Figure 2). 
'These two frequency channels have been combined with- 
out taking the differences in passbands into account, 
leading to CO(1-0) being emphasized on large Planck- 
dominated scales in the coadded map, but not on small 
ACT-dominated scales. This likely explains the haziness 
of the emission in purple in Figure 4, where the low- 
frequency channel (red) contains significant CO(1-0) 
emission in the Planck map but is dominated by other, 
less prominent emission mechanisms in the ACT map. 
A quantitative interpretation of the frequency spectra of 
particular regions in the Galactic center will therefore re- 
quire careful consideration of bandpass effects, and pos- 
sibly the use of external spectroscopic data (e.g., Dame 
et al. 2001; Eden et al. 2020) and/or the CO component 
maps from Planck (Planck Collaboration X 2016). Such 
spectral analysis will be the subject of future work, and 
for now we urge caution when interpreting the colors in 
Figure 4 in terms of emission mechanisms or spectral 
indices. 


5. POLARIZATION MAPS 


Figure 3 presents the full-resolution Stokes Q and U 
maps obtained through the mapmaking algorithm at 
each frequency band. A striking feature of the maps 
is the strong polarization signal of the GCRA, extend- 
ing roughly from b — —0.5? to b — 0.5? in both f090 and 
f150. The signal is weaker in f220, which is dominated by 
polarized dust emission. Strong polarized signals can be 
generally seen near the CMZ along the Galactic plane 
across all frequency bands, with especially prominent 
polarization features near regions such as Sgr A* and 
Sgr B2. This suggests that the observed polarization 
signals are not dominated by diffuse emission along the 
LOS, but rather by emission directly from the CMZ. 
Since we are focusing on high S/N regions (= 3) that 
are negligibly impacted by debiasing, we do not debias 
the polarization quantities (Plaszczynski et al. 2014). 

'To create a three-color polarization image analogous 
to that in total intensity, we first compute the polarized 
intensity P = yQ? + U? in each band. We synthesize 
the three polarized intensity maps into a three-color im- 
age using f090, f150, and f220 as the red, green, and 
blue channels, respectively. The result is shown in the 
lower panel of Figure 4. The polarized emission has a 
strikingly different morphology than total intensity (see 
upper panel of Figure 4). The polarized GCRA stands 
out distinctively from the background in red, indicat- 
ing dominance of f090, consistent with the prominence 
of synchrotron radiation in this region. Similarly, ra- 
dio sources including Sgr A* and G359.23-0.82 ("the 
Mouse"; e.g., Predehl & Kulkarni 1995) appear red. On 
the other hand, molecular cloud complexes such as the 
l = 1.3 complex (Bally et al. 1988) and G0.55-0.85 
(RCW 142; Gardner & Whiteoak 1975) appear blue, 
consistent with the expected predominance of dust emis- 
sion. 

One quantity of interest is the polarization angle, de- 


fined as 
1 U 
p = 5 arctan (5) j (4) 


'The polarization angle is directly related to the plane- 
of-sky magnetic field orientation by a 90? rotation. Dust 
grains tend to align their short axes parallel to the mag- 
netic field, while they radiate photons preferentially po- 
larized parallel to their long axes. The synchrotron 
polarization angle, or electric vector position angle, is 
similarly orthogonal to the local magnetic field orienta- 
tion for optically thin emission. Hence, the magnetic 
field orientation is orthogonal to the polarization angle 
in both emission mechanisms. We note, however, that 
dust and synchrotron emission do not necessarily trace 
the same magnetic field, as they generally probe differ- 
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Figure 6. A visualization of magnetic field orientations using line-integral-convolution (LIC) with a 1? kernel. Contours in 
the map trace magnetic orientations. Rows represent f090, f150, and f220 respectively. Total intensity maps are shown in the 


background with the same color scales in Figure 1. 


ent volumes along the LOS. The observed magnetic field 
morphology at a given frequency depends on the relative 
contribution of different emission components, which in 
turn depends on the spatial distribution of dust density 
versus cosmic ray density and the underlying magnetic 
field orientation and strength (see Han (2017) for a re- 
view). 

Figure 6 presents a visualization of the inferred mag- 
netic field orientation in each of our bands using line 
integral convolution (LIC; Cabral & Leedom 1993) with 
a kernel size of 0.5?. Each contour in the map traces 
the magnetic field orientation. The magnetic field is ap- 


proximately parallel to the Galactic plane near the CMZ 
for both f090 and f150, and is noticeably tilted for f220 
within the range |I| S 1.5°. In particular, within a box 
of |I| < 1.5?, |b] < 0.15? we measure the mean polar- 
ization angle to have a tilt of = 20? with respect to the 
Galactic plane, consistent with the ~ 22? tilt previously 
noted by, e.g., PILOT (Mangilli et al. 2019). 

'The f090 map is noticeably more disordered, with es- 
pecially prominent features at the GCRA, where the 
plane-of-sky magnetic field is aligned with the orienta- 
tion of the arc. This 90? flip in polarization angle at the 
GCRA has been observed by the QUIET Collaboration 
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Figure 7. Polarization fractions (background) and magnetic field orientation (line segments) are shown for our three bands 
(f090, f150, and f220). To estimate the magnetic field orientations, the polarization field is smoothed with a Gaussian kernel 
FWHM=2", and then resampled with a pixel size of 2'. Line segments with large uncertainty in polarization angle 67 > 15° are 
masked. 
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(Ruud et al. 2015) at both 43 GHz and 97 GHz. This or- 
thogonal feature is less prominent at f150 and disappears 
at f220, as expected from a synchrotron-dominated sig- 
nal. 

The polarization fraction p = yQ? + U?/I in each 
band is shown in Figure 7. In each panel, we overlay 
the magnetic field orientation in the CMZ at 2' resolu- 
tion. Along the Galactic plane the polarization fraction 
is generally low, p S 296. This is consistent with the 
previous observations from, e.g., Planck (Planck Col- 
laboration et al. 2020b) and PILOT (Mangilli et al. 
2019) that found polarization fractions at the percent 
level (< 1.5%) in the Galactic center region. We see co- 
herent magnetic fields even within regions of relatively 
low polarization fraction, in agreement with both cloud- 
scale observations and the relatively few wide-area dust 
polarization measurements, both of which tend to find 
very ordered magnetic fields (Chuss et al. 2003; Pillai 
et al. 2015; Mangilli et al. 2019). The large-scale co- 
herence in the inferred magnetic field direction suggests 
that the polarized emission is dominated by the CMZ. 
The low polarization fraction could be due to one of sev- 
eral effects, or to a combination of them. Perhaps the 
most likely is that the magnetic field orientation fluctu- 
ates both along the LOS and within the beam smoothing 
radius, resulting in depolarization. There are so many 
emitting regions along the LOS in the Galactic disk that 
small variations in the magnetic field orientation aver- 
age out in the LOS integration, such that observed de- 
viations from the mean magnetic field orientation are 
small. We note, however, that simulations of the Galac- 
tic magnetic field used to interpret PILOT data sug- 
gest that this effect may not be sufficient on its own to 
account for the entirety of the observed depolarization 
(Mangilli et al. 2019). Another possibility is that the 
mean field has a significant LOS component. Because 
magnetically aligned dust grains spin around their short 
axes, the net dust emission is more strongly polarized for 
regions with a predominantly plane-of-sky magnetic field 
than for regions where the magnetic field is more par- 
allel to the LOS. However, a significant LOS magnetic 
field component would not be expected to dominate the 
entirety of the CMZ if the magnetic field has a signifi- 
cant azimuthal component. Finally, it may be that the 
mean field in the CMZ is itself a product of superim- 
posed, misaligned structures that each have large-scale 
coherence, e.g., the twisted ring geometry proposed for 
the distribution of dust density in the CMZ (Molinari 
et al. 2011). While possible, such a scenario demands 
great uniformity in the relative total and polarized in- 
tensities in each component to avoid dispersion in the 
observed polarization angles. On balance, we favor a 


coherent magnetic field in the CMZ dust, with LOS dis- 
order as the primary driver of low polarization fractions, 
but more detailed modeling of the present data is war- 
ranted to assess the relative importance of each of these 
effects. 


6. NOTABLE OBJECTS 


With arcminute resolution in three frequency bands, 
we detect many known radio and infrared sources, some 
of which have not been previously observed at ACT 
frequencies. Although the main focus of this paper is 
presentation of the Galactic center coadded maps, in 
this section we demonstrate the fidelity of these maps 
and their broad potential for different scientific investi- 
gations by highlighting select objects. All objects dis- 
cussed in this section are marked in Figure 5, which in- 
cludes additional selected radio sources listed in LaRosa 
et al. (2000) and submillimeter sources from the CM- 
Zoom Survey (Battersby et al. 2020) visible in our maps. 
This list of notable sources is non-exhaustive, and in 
particular, our maps extend to a wider range in Galac- 
tic longitude than either the LaRosa et al. (2000) or 
Battersby et al. (2020) catalogs. 


6.1. Sgr A and GCRA 


Sagittarius A (Sgr A) is a complex radio source located 
at the center of our Galaxy. It consists of Sgr A East, an 
extended nonthermal source with a radius of — 3', and 
a thermal source Sgr A West, which has three-arm spi- 
ral morphology and lies within Sgr A East (e.g., Ekers 
et al. 1983; Yusef-Zadeh & Morris 1987b; Anantharama- 
iah et al. 1991). Infrared monitoring of stellar orbits in 
the vicinity of Sgr A has also revealed the existence of 
a supermassive black hole Sgr A* that lies within Sgr A 
West (e.g., Ghez et al. 2008) and acts as the dynamical 
center of our Galaxy (Backer & Sramek 1999). 

The region of sky surrounding Sgr A* has been the 
subject of extensive multifrequency observations both 
in imaging and polarimetry (e.g., Stolovy et al. 1996; 
Bower & Backer 1998; Melia et al. 2000; Baganoff et al. 
2003; Chuss et al. 2003). Polarized observations in the 
millimeter bands, in particular, are important for under- 
standing the accretion process near the black hole and 
associated relativistic emission (e.g., Agol 2000; Melia 
et al. 2001). Linear polarization of Sgr A* at millime- 
ter wavelengths was first reported by the Submillimetre 
Common-User Bolometer Array (SCUBA; Aitken et al. 
2000), which they interpret as synchrotron-dominated 
polarized emission sourced by the gas in the vicinity of 
the black hole. The observed polarization fraction of 
Ser A* is ~ 3% at 2 mm. Subsequent interferomet- 
ric imaging surveys (e.g., Macquart et al. 2006; Mar- 
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Figure 8. GCRA and Sgr A*. The left panel shows the polarized intensity in the region, measured from f090 coadded. Contours 
show levels of total intensity at f090 with a spacing of 2 MJysr”* up to 30 MJysr”*. The right panel shows the inferred magnetic 
field orientations from the f090 map as line segments in 0.5' pixelization (full resolution). Segments are shown with varying 
opacity that scales linearly with the S/N in polarized intensity and saturates when S/N = 3. In the background we show a 
radio image of the region from MeerKAT (Heywood et al. 2019) which observes at 1.28 GHz in 6" pixelization. The expected 
location of Sgr A* is indicated with a white cross mark in both panels. Note that the MeerK AT image is shown for visualization 
purposes only, as no primary beam corrections have been applied, and the entire Galactic plane is seen through the primary 
beam sidelobes. Caution should be taken when interpreting the numerical values in this image (see Heywood et al. (2019) for a 


detailed discussion). 


rone et al. 2006) measured a ~ 2% polarization frac- 
tion at 3.5 mm, and larger values at higher frequen- 
cies. Strong emission centered on Sgr A* is visible in 
the coadded maps, showing up clearly in the multifre- 
quency image with a yellow color in total intensity (see 
the upper panel in Figure 5), implying a predominance 
of synchrotron emission in the region. Its location in- 
dicates that the emission is likely dominated by Sgr A* 
itself instead of the overlapping components in Sgr A 
that are unresolved with the ACT beam. Regions sur- 
rounding Sgr A* are polarized at 2 — 496 level, as seen 
in Figure 7 for f090 and f150, and show up as a reddish 
“blob” in the multifrequency polarimetry (see the lower 
panel in Figure 5). This may be due to synchrotron 
emission from the nearby nonthermal filaments within 
a beam smoothing radius. The polarized emission in 


the vicinity of Sgr A* has a lower polarization fraction 
of ~ 1.5% at all three bands, consistent with the de- 
polarization noted by SCUBA (Aitken et al. 2000) at 
2 mm. The slightly lower polarization fraction seen in 
the ACT data is likely due to a beam depolarization ef- 
fect from the larger ACT beam (^ 2”) in comparison to 
the SCUBA beam (^ 34" at 2 mm). 

In Figure 8 we present a zoom-in view of the region 
surrounding Sgr A*. The left panel shows the polarized 
signal in f090 overlaid with contours from the total in- 
tensity in f090. Strong emission from Sgr A* is seen in 
total intensity but not in polarization, where the emis- 
sion is more diffuse and extends ~ 3' away from the 
central source. This is further evidence that the polar- 
ized signal in the vicinity of Sgr A* is emitted by the 
surrounding non-thermal filaments, while the emission 
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Figure 9. Molecular cloud known as “the Brick”. Left: total intensity measured from ACT+Planck f220 coadd map is plotted 
in the background. The Herschel 500 ym measurements (Molinari et al. 2016) are shown as contours indicating 50th, 70th, 
and 90th percentiles from lighter to darker contours. Right: total intensity measured by Herschel 500 ym is shown in the 
background. We show the magnetic field orientation inferred from the f220 map as line segments. Segments are shown with 
varying opacity that scales linearly with the S/N in polarized intensity and saturates when S/N = 3. 


from Sgr A* itself is highly depolarized. In the right 
panel we show the inferred magnetic field orientations 
from the polarized signal at f090 overlaid on top of a 
radio image of the same region from MeerKAT (Hey- 
wood et al. 2019), which observes at 1.28 GHz with a 
6” beam. The magnetic field morphology inferred from 
our f090 map closely follows the underlying non-thermal 
filamentary structure. The morphology is also in broad 
agreement with previous Caltech Submillimeter Obser- 
vatory (CSO; Chuss et al. 2003) observations at a wave- 
length of 350 um with a 20” beam. 

Figure 8 also shows the GCRA, a prominent radio fea- 
ture located at ~ l = 0°10’, which consists of a bundle 
of thin filaments running perpendicular to the Galactic 
plane (e.g., Yusef-Zadeh & Morris 1987a; Ananthara- 
maiah et al. 1991). The GCRA is known to be a highly 
polarized synchrotron source, though its origin is still 
poorly understood. The strong synchrotron emission 
implies that free electrons are present in the GCRA and 
are accelerated to relativistic speeds in the presence of 
a strong magnetic field in the region. Various models 
have been proposed to explain the source of electrons 


and the acceleration mechanism (see, e.g., Serabyn & 
Morris 1994, for a review), though the matter is still 
under debate. 

In millimeter bands, the GCRA has previously been 
detected at 7 mm (Reich et al. 2000), 3 mm (Pound & 
Yusef-Zadeh 2018), and 2 mm (Staguhn et al. 2019), 
which the latter notes was the highest-frequency de- 
tection of the GCRA at the time. Polarized emission 
from the GCRA has also been previously detected at 2 
and 3 mm by Culverhouse et al. (2011), and at 3 mm 
and 7 mm by Ruud et al. (2015). In our coadded 
maps, GCRA appears in total intensity in both f090 
and f150. The associated polarized emission can also be 
seen clearly in f090 and f150 with polarization fractions 
reaching ~ 30%. This is considerably higher than the 
~ 10% peak polarization noted by the QUaD Galactic 
Plane Survey (Culverhouse et al. 2011) at the same fre- 
quencies, likely due to the improved angular resolution 
in our coadded maps (2' at f090, 1.4’ at f150) in com- 
parison to Culverhouse et al. (2011) (5' at 100 GHz, 3.5' 
at 150 GHz). The polarized emission from the south- 
ern portion of the GCRA is also visible in f220, which 
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Figure 10. A cloud triad known as “the Three Little Pigs” consisting of G0.145-0.086 (“Straw Cloud”), G0.106-0.082 


(“Sticks Cloud"), and G0.068-0.075 (“Stone Cloud"). 


The data are plotted following Figure 9, with the left panel showing 


the ACT+Planck f220 map with the Herschel 500 um image overlaid as contours (indicating 50th, 70th, and 90th percentiles 
from lighter to darker colors), and the right panel showing the Herschel 500 ym map with the magnetic field orientations inferred 
from the f220 map overlaid as line segments. Segments are shown with varying opacity that scales linearly with the S/N in 


polarized intensity and saturates when S/N= 3. 


is likely the highest frequency at which this structure 
is detected to date (note especially the {220 Q map in 
Figure 3). In addition to being fainter at 220 GHz on 
account of the falling synchrotron spectrum, the GCRA 
is also obscured by emission from dust along the LOS. 
The uniformity of the polarized emission observed in 
the Arc as seen in Figure 3 implies that a highly or- 
dered magnetic field exists along the Arc that deviates 
sharply from the large-scale magnetic field geometry (see 
Figure 6). In particular, the magnetic field orientation 
inferred from f090 (as seen in the right panel of Figure 8) 
aligns closely with the filamentary structure perpendic- 
ular to the Galactic plane. This is in broad agreement 
with the morphology observed at 43 GHz and 96 GHz 
by QUIET with lower angular resolution (Ruud et al. 
2015). 


6.2. The Brick 


G0.253+0.016, also known as “the Brick”, is a dense, 
massive molecular cloud in the CMZ, and a prominent 
infrared dark cloud (Carey et al. 1998; Longmore et al. 
2012). In the context of understanding the low star 
formation rate in the Galactic center environment, the 
Brick is a particularly interesting case study. Despite its 


high mass (> 107 M¿) and density (> 10* cm-?), evi- 
dence of star formation is nearly absent in the Brick, and 
thus it may provide an ideal opportunity to study the 
initial conditions of high-mass star formation (Lis et al. 
1994; Longmore et al. 2012; Kauffmann et al. 2013; Mills 
et al. 2015; Walker et al. 2021). A number of factors have 
been invoked to explain the dearth of star formation in 
G0.2534-0.016, including solenoidal turbulence driven by 
strong shear in the CMZ (Federrath et al. 2016; Krui- 
jssen et al. 2019; Dale et al. 2019; Henshaw et al. 2019), 
or strong cloud-scale magnetic fields (B ~ mG, Pillai 
et al. 2015). 

The Brick stands out at high contrast to the back- 
ground in the coadded total intensity maps at both 
150 and 220GHz. Our polarization measurements at 
these frequencies probe the magnetic field structure in 
the dust toward G0.253+0.016 at ~arcminute scales. 
These observations complement 20" resolution polariza- 
tion data at 350 um from the Caltech Submilimeter Ob- 
servatory (CSO; Dotson et al. 2010; Pillai et al. 2015). 
We find that the inferred magnetic field orientation is 
aligned parallel to the long axis of the Brick on the plane 
of the sky (Figure 9), and the polarization angles are 
very ordered in this region, in agreement with the CSO 
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Figure 11. G359.23-0.82 or “the Mouse” is a pulsar wind nebula (PWN) traveling with high velocity (7:300 kms”*) with 
respect to ISM, causing a comet-like tail. The left panel shows the total intensity in f090 with magnetic field orientation 
over-plotted in line segments. Both the background and magnetic field are smoothed to a resolution of 2.2” to increase the 
signal-to-noise ratio. Segments are shown with varying opacity that scales linearly with the S/N in polarized intensity and 
saturates when S/N= 3. The middle panel shows the polarized intensity in f090 after smoothed to a resolution of 2.2”. The right 
panel shows a radio image of the region from MeerKAT (Heywood et al. 2019) which observes at 1.28 GHz in 6” pixelization, 
with the magnetic field orientation from f090 over-plotted as line segments similar to the leftmost panel. 


data at smaller angular scales. Pillai et al. (2015) use the 
strong coherence of the magnetic field orientation in the 
Brick to compare the inferred magnetic field strength to 
the gas velocity dispersion measured from N2H* emis- 
sion (Kauffmann et al. 2013). Those authors find that 
magnetic fields dominate over turbulence in the Brick. 
The coherent magnetic field structure in our observa- 
tions is consistent with the expectation that turbulence 
in the Brick is sub-Alfvénic at the scales probed by ACT. 
The ACT polarized emission is brightest at the northern 
part of the Brick, with a peak polarization fraction of 
1.896. The polarized intensity is lower in the southern 
portion of the cloud, and the SNR on the polarized in- 
tensity drops below 3. This depolarization may be due 
in part to unresolved polarization structure within the 
ACT beam, and/or to incoherent contributions to the 
polarized emission along the LOS. 


6.3. The Three Little Pigs 


The cloud triad G0.145-0.086, G0.106-0.082, and 
G0.068-0.075 visible in Figure 10 has been dubbed “the 
Three Little Pigs. All three clouds have been noted 
as a set of compact dusty sources in the CMZoom Sur- 


vey (Battersby et al. 2020), while G0.068-0.075 also ap- 
pears in the SCUBA-2 Compact Source Catalog (Par- 
sons et al. 2018). As Figure 10 illustrates, each cloud is 
also apparent in the 500 jum data from Herschel Infrared 
Galactic Plane Survey Herschel (Hi-GAL; Molinari et al. 
2016). Interestingly, the 3" resolution 230 GHz obser- 
vations with the Submillimeter Array as part of the 
CMZoom Survey have revealed a dearth of substruc- 
ture in G0.145-0.086 (“Straw Cloud"), somewhat more 
substructure in G0.106-0.082 (“Sticks Cloud"), and yet 
more in G0.068-0.075 (*Stone Cloud"). 

ACT f220 measurements give a first look at the mag- 
netic field geometry in these clouds at arcminute res- 
olution. The Straw Cloud, perhaps owing to a lower 
column density or lack of dense substructure, has a mag- 
netic field orientation that deviates little from the large- 
scale field structure. In contrast, both the Sticks and 
Stone Clouds have polarization angles in their interiors 
that are highly misaligned with the large-scale magnetic 
field. Similar to the depolarization observed toward the 
Brick, the cancellation of polarized emission from dust 
in different regions within the cloud and/or other dust 
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Figure 12. G357.7-0.1, or “the Tornado”, is typically classified as a supernova remnant. The left plot shows the total intensity 
in its neighborhood in f090 coadded map. Line segments indicate The right plot shows the magnetic field orientations inferred 
from f090. They are shown with varying opacity that scales linearly with the S/N in polarized intensity and saturates when 
S/N= 3. The right panel shows the corresponding polarized intensity map in f090. Both maps are shown at the full resolution 


from mapmaking (0.5’). 


along the LOS may explain the low polarized intensities 
observed, particularly in the Stone Cloud. 


6.4. The Mouse 


G359.23-0.82, also known as “the Mouse”, is a PWN 
powered by the young X-ray source PSR J1747-2958 
(Predehl & Kulkarni 1995; Camilo et al. 2002). 
G359.23-0.82 was originally discovered in radio contin- 
uum data from the Very Large Array (VLA), and derives 
its nickname from its bright compact nebula “head” 
and extended radio “tail” (Yusef-Zadeh & Bally 1987). 
The Mouse is strongly linearly polarized at centime- 
ter wavelengths (Yusef-Zadeh & Gaensler 2005). Dis- 
tances to PSR J1747-2958 and the Mouse are uncer- 
tain, but they are not at the Galactic center: observa- 
tions of neutral hydrogen absorption set the maximum 
distance to G359.23-0.82 at ~ 5.5 kpc (Uchida et al. 
1992). Gaensler et al. (2004) argue for a distance of 
~ 5kpc, a value now commonly adopted (e.g., Klin- 
gler et al. 2018). At 5kpc, the transverse velocity of 
PSR. J1747-2958 is 306 + 43 kms"! (Hales et al. 2009). 
'The Mouse is a striking example of a bow shock nebula, 


formed by the interaction of the pulsar with the ambient 
ISM as it travels at supersonic speeds (e.g., Gaensler & 
Slane 2006). 

The Mouse is a prominent object in the ACT f090 
map, both in total and polarized intensity (Figure 11). 
In particular, polarized emission is detected significantly 
across the peak of the Mouse, which is expected for 
a PWN. Significant polarized emission is also detected 
along its tail, and exhibits a similar morphology as seen 
by MeerKAT at 1.28 GHz (Heywood et al. 2019) with 
a 6" beam, albeit at lower resolution in the ACT data. 
The implied magnetic field orientation in the f090 band 
is roughly parallel to the Mouse's extended tail, con- 
sistent with observations at 3.5 and 6 cm by the VLA 
(Yusef-Zadeh & Gaensler 2005). The Mouse is travel- 
ing eastward in declination, which is roughly toward the 
lower left-hand corner of Figure 11. 


6.5. The Tornado 


G357.7-0.1, “the Tornado,” is typically classified as 
a SNR, though its unusual properties have prevented 
a definitive explanation (Gaensler et al. 2003; Chawner 
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Figure 13. l = 1.3 molecular complex. The left plot shows the total intensity in f220 (smoothed with FWHM=1’) with contours 
indicating the 50th, 70th, 90th percentiles in the Herschel 500 ym map. The right plot shows the Herschel 500 ym map with 
magnetic field orientation inferred from the f220 map as an overlay, after smoothed to a resolution of 1.4’. Segments are shown 
with varying opacity that scales linearly with the S/N in polarized intensity and saturates when S/N= 3. 


et al. 2020). The Tornado has been long observed in 
radio imaging and polarimetry (e.g., Milne 1970; Shaver 
et al. 1985; Law et al. 2008), which consistently show 
a bright “head” region and a “tail” region roughly 10" 
in extent. Recently, mid- and far-infrared dust emission 
has been detected with Spitzer and Herschel, revealing 
a large dust reservoir in the head region ( 17 Mọ) and 
consistent with interstellar matter swept up in a super- 
nova blast wave (Chawner et al. 2020). The head of the 
Tornado has also been detected by Chandra in X-rays 
without evidence for embedded point sources (Gaensler 
et al. 2003), lending further credence to its classifica- 


tion as a SNR. However, the provenance of the tail is 
still unresolved (see Chawner et al. 2020, for a recent 
discussion). 

'The Tornado is prominent in the f090 and f150 Stokes 
Q and U maps, but not f220 (see Figure 3). Likewise, 
the region stands out in reddish brown in the three-color 
polarization map (Figure 5). This suggests the promi- 
nence of synchrotron emission in this source. A closer 
examination of the Tornado in the f090 band is presented 
in Figure 12. Here, we see the extended tail region in 
total intensity but not in polarization, while the head 
is prominent in both. This morphology is consistent 


with 4.9 GHz polarimetric observations by Shaver et al. 
(1985). The inferred magnetic field at f090 is approxi- 
mately perpendicular to the extended tail in the eastern 
side of the Tornado, and is tilted toward the head on 
the western side. This is also in a broad agreement with 
the magnetic field morphology noted by (Shaver et al. 
1985) at 4.9 GHz. We observe a maximum polarization 
fraction of the Tornado in f090 of 8.5% + 1%, slightly 
lower than the ~ 10% observed at 4.9 GHz at signifi- 
cantly higher resolution (12 x 26" beam, Shaver et al. 
1985). It is likely that much of the difference is due to 
more beam depolarization in the ACT data. 


6.6. 1 = 1.3 Complex 


The combination of ACT and Planck data used in the 
coadded maps enables large regions to be mapped with 
fidelity on both large and small angular scales. Likewise, 
the high sensitivity of the polarimetry permits mapping 
of more diffuse regions of molecular clouds, not just 
bright cores. These capabilities are highlighted in the 
20 x 30’ maps of the | = 1.3 complex in Figure 13. 

The | = 1.3 complex is a large, high-velocity- 
dispersion molecular cloud complex extending from 
roughly 1.2-1.6? in Galactic longitude (Bally et al. 
1988). The elevated abundance of SiO and high ratio 
of CO(3-2) to CO(1-0) emission in some clouds within 
the complex suggest the presence of strong shocks, per- 
haps from cloud-cloud collisions or supernova explosions 
(Huettemeister et al. 1998; Oka et al. 2001; Rodriguez- 
Fernandez et al. 2006; Tanaka et al. 2007; Parsons et al. 
2018; Tsujimoto et al. 2021). This complex may sit at 
the intersection of a dust lane with the nuclear ring, sup- 
plying it with material (Huettemeister et al. 1998; Fux 
1999; Rodriguez-Fernandez et al. 2006; Liszt 2008). 

Total emission from the l = 1.3 in f220 and Herschel 
500 um (Molinari et al. 2016) is presented in Figure 13, 
with good morphological correspondence between the 
two maps. In the right panel, we overlay the f220 mag- 
netic field orientation on the higher resolution Herschel 
map. While many density structures show clear align- 
ment with the magnetic field orientation, this is not uni- 
versally observed. The highest intensity regions have 
comparatively low polarized intensities, suggesting ele- 
vated magnetic field disorder or a loss of grain alignment 
in the densest regions. 


7. SUMMARY AND FUTURE PROSPECTS 


We have presented new arcminute-resolution maps of 
the Galactic center region at microwave frequencies by 
combining data from ACT and Planck. Known radio 
features appear at high significance in both total inten- 
sity and polarization in three frequency bands. The po- 
larization maps provide a frequency-dependent probe of 
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magnetic fields, demonstrating a change in the observed 
magnetic field morphology as the fractional contribu- 
tions of synchrotron radiation and thermal dust emission 
from different regions within the Galactic center along 
the LOS vary with frequency. With wide-field maps 
at higher angular resolution, we identified known radio 
sources and molecular clouds, some of which have not 
previously been observed in polarization at microwave 
frequencies. With three frequency bands, our total in- 
tensity maps reveal the rich physical environment in the 
CMZ with spatially varying combinations of different 
emission mechanisms, including synchrotron, free-free, 
dust, and molecular line emission in the CMZ. Separa- 
tion of these emission components will be the subject of 
a follow-up paper. 

The coadded maps produced in this work are now 
publicly available on the NASA Legacy Archive Mi- 
crowave Background Data Analysis (LAMBDA; Miller 
& LAMBDA 2018)?. These maps are suitable for trac- 
ing magnetic field morphology across the Galactic center 
region and measuring the total and polarized emission 
from individual sources. However, caution is urged for 
multifrequency analyses due to the bandpass mismatch 
between ACT and Planck that results in a slight scale- 
dependence of effective band centers for different emis- 
sion mechanisms. As discussed in Section 4, CO(1-0) 
emission falls within the Planck 100 GHz passband but 
not f090, amplifying bandpass mismatch effects in the 
resulting coadded map. 

ACT has continued to observe the Galactic center dur- 
ing 2020, collecting a similar amount of data to that 
used in this work. In addition, the daytime data from 
both 2019 and 2020 can, in principle, be corrected for 
thermal telescope distortions (Aiola et al. 2020), which 
would again double the total amount of data. There- 
fore, ACT maps with half the pixel noise variance of 
those presented here are possible based solely on data 
that has already been collected. Additionally, we plan 
to apply the mapping techniques used here to approxi- 
mately 70 degrees of the Galactic plane covered by ACT 
from 2019 and 2020. Furthermore, the addition of the 
low-frequency array to ACT in 2020 (Li et al. 2018; Si- 
mon et al. 2018) will also allow us to map the Galactic 
plane at 27 GHz and 39 GHz, likely yielding new insights 
on the Galactic center environment. 

'The next observational step at these frequencies will 
be the Large Aperture Telescope of the Simons Obser- 
vatory (Ade et al. 2019), anticipated to see first light in 
2023 from the same site in Chile. This new instrument 


2 https: / /lambda.gsfc.nasa.gov/product /act/actadv. sr. gc. 1. info. 
cfm 
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will have the same 6-meter diameter primary as ACT, 
but with an instrumented focal plane of 5 times larger 
area (Zhu et al. 2021). The nominal scan strategy will 
continuously cover the entire sky in the declination range 
between 4-25? and —40°, providing coverage of over 100 
degrees of the Galactic plane in five frequency bands 
in both total intensity and polarization. The five-year 
map noise should improve on ACT by roughly a factor of 
three. The Galactic center will be observed at higher fre- 
quencies by the CCAT-prime project (Choi et al. 2020) 
and will also be a good target for future balloon-borne 
instruments, which can achieve sub-arcminute resolution 
with similar sensitivity at even higher frequencies, e.g., 
BLAST Observatory (Lowe et al. 2020). By 2030, we 
can also anticipate data from CMB-S4 (CMB-S4 Col- 
laboration et al. 2016), with an additional map noise 
improvement by a factor of four. This unrivaled com- 
bination of resolution, sky coverage, and sensitivity at 
microwave frequencies will enable many new inquiries 
into the properties of the Milky Way. 
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APPENDIX 


A. CALIBRATION METHOD 
In ACT DR4 (Aiola et al. 2020), the data were calibrated from raw data acquisition units to physical units with 


dPW = qP^9 x Ras x gatm X f, (A1) 


where dP“, ¿PAC represent detector data in physical unit and data acquisition unit respectively, Rgs represents the 
intrinsic responsivity of each detector measured from the most recent bias step, Jatm is the atmospheric correction 
factor and f represents an optical flatfield. Both gatm and f are estimated from detector responsiveness to the 
atmospheric signal which is treated as a common mode for all detectors. In the presence of non-atmospheric thermal 
contamination signatures this approach to estimating the common mode may bias the calibration. Preliminary analyses 
on the Advanced ACT data collected after 2017 have shown evidences for the presence of such thermal contamination. 
Hence, the calibration method needs to be updated to account for this bias. 
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To circumvent the problem, we switch from a common-mode based calibration to a planet based approach, given by 


dPW = qP^9 x Rpg x fo, (A2) 


where we have dropped the atmospheric correction gat. fp now represents an optical flatfield measured from detector 
responsiveness to emission from Uranus instead of the atmosphere. This leads to an improved calibration model and 
better gain stability, and this is expected to be the standard calibration method for future data releases of ACT. 

One caveat to note is that ACT DR4 (and earlier) maps are calibrated to Planck maps in a final step. This step, 
however, is missing here due to the preliminary nature of the Advanced ACT data in 2019, but since the Planck 
calibration is found to be close to the planet calibrations for ACT DR4, and we are using an improved calibration 
model from DRA, we do not expect this to be a concern and estimate a O(1%) uncertainty in global gain calibration 
as a result. 


B. BEAM LEAKAGE CORRECTION 


'To reduce the contamination in the polarization maps caused by beam leakage effects, we modeled the observed 
polarization maps Pops as a sum of a beam-convolved sky map Psky and a leakage component, given by 


Pobs(x) = Pay (x) + I on Toy (k) Bog (k)e"™**, (B3) 


where the leakage component is a convolution of the beam-convolved temperature map Tsxy with an effective beam 
Beg given in Fourier space by 
Ber (k) = Bp*(k)Bp(k) . (B4) 


This represents a combination of deconvolving the temperature beam Br and convolving with a leakage beam Bp. 
Our strategy is then to build a model of Bes(k), convolve it with the temperature map, and deproject it from the 
observed polarization maps for each dataset. 

As we assumed Uranus signal is unpolarized, any signal we measure in the polarization is a sign of T-to-P leakage. 
Hence, we modeled Bog using observations of Uranus made in the same observation season, with the following steps: 


1. We made Uranus planet maps for each Uranus observation in a source-centered reference frame with the scan 
direction as the horizontal axis. 


2. We reprojected each planet map into the Galactic coordinate system. As the Q/U reference frame needs to be 
rotated depending on the scan directions, and the Galactic center observations consist of two scan directions 
taken during rising and setting respectively, we accounted for the difference in scan directions by rotating the 
Q/U reference frame for both rising and setting scan directions and performed a weighted average (for each planet 
map) depending on the total number of rising and setting scans made during Galactic center observations. 


3. We stacked all reprojected Uranus maps for each dataset (per frequency band per array) with inverse-variance 
weighting to obtain estimates of Br and Bp for each dataset. 


4. We performed a real space cut to remove noise outside a radius of rmax in Br and Bp for each dataset, then 
calculated Beg (k) = Bp/Br in Fourier space. 


5. We further cleaned Bog (k) with a k-space cut k < kmax to get rid of small scale noise in the beam model. We 
then re-filled the k-space outside kmax by mirroring the value at k = kmax with Beg (k)|k>kmax = Bog (k kmax/ k): 
In practice, the details of the extrapolating function make little difference. This specific function is chosen to 
ensure that the transition at k = kmax is smooth, and it extends naturally to infinity. 


In Table 2, we listed the choices of rma, and kmax for each dataset. As a result of these steps, we obtained Bog(k) for 
each dataset respectively. Treating the coadded temperature map as the “true” sky model T3, at each frequency band, 
we predicted the expected T-to-P leakage from the derived leakage beam model Beg (k) and subsequently deprojected 
the expected leakage from each observed map. The resulting ACT maps after leakage correction are shown in Figure 14 
in comparison to the coadded maps. 
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Figure 14. Comparison between ACT-only maps (left column) and ACT+Planck coadded maps (right column) in total 
intensity, similar to Figure 1. Rows from top to bottom correspond to f090, f150, and f220, respectively. The ACT maps shown 
here are produced by coadding maps from different detector arrays with inverse variance weightings at each frequency band 
respectively. Each map is plotted on a logarithmic color scale from 0.3-30 MJysr”* for f090 and f150, and from from 3-100 
MJysr”* for £220. 


array freq Tmax kmax 
PA4 f150 5 14500 
PA4 f220 4' 20000 
PA5  f090 g 10500 
PA5 f150 5 16500 
PA6  f090 g 10500 
PAG  f150 5 16500 


Table 2. Parameters used when building a 2D leakage beam model (in step 4 and step 5). 
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C. COADDED F220 MAPS 


Maps at f220 are much noisier than the two other bands, so using the same coadding pipeline results in a slow 
convergence. Instead, we adopted a simpler Fourier-based coadding algorithm. Denote the ACT f220 map as mı, and 
Planck 217 GHz map as mg. The coadd map Meoaaa was obtained with a simple inverse-variance coadding as 


Meoadd = (Ni + N 5D” (Nim + Nz m3) P (C5) 


where Nj and Nə are noise covariance matrices assumed for ACT and Planck maps respectively. Specifically, we 
assumed simple Fourier-based noise model with Ny = w?{1 + (€/€cnee)~*] and Na = w2b;/b», with w, = 268.97 uK’, 
W = 124.22 uK' the noise levels of ACT and Planck maps respectively, lknee = 4000, and a = —3. bı and ba represent 
the beam model in ACT and Planck respectively, and the factor bı /b2 in N2 represents the effect of a combination of 
deconvolution of the Planck beam and convolution with the ACT beam in the Planck noise model. In addition, we 
further applied a high-pass Butterworth filter in the ACT map with Ze = 2500 and a = —5 along the two cross-linked 
scan directions with a width given by the beam FWHM in Fourier space. This helps suppress excess noise along the 
scan directions. We then applied an additional high-pass Butterworth filter with £e = 200 and a = —10 to suppress 
the large-scale atmospheric noise in the ACT f220 map. Both of these filters are included in the ACT noise model N;. 

'The inverse-variance map of the output coadd map was estimated using the same weighted average from Planck 
inverse-variance map V, ! and ACT inverse-variance map M3 !, given by 


Neoadd = GNM, + CN (C6) 


with ci, c» defined as the -space mean of c1, = Nia (NL! + NT from £ — 2500 to 5000 where we expect to be 
dominated by white noise. 
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